Decoherence of a quantum gyroscope 
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We study the behavior of a quantum gyroscope, that is, a quantum system which singles out a 
direction in space in order to measure certain properties of incoming particles such as the orienta- 
tion of their spins. We show that repeated Heisenberg interactions of the gyroscope with several 
incoming spin- 1/2 particles provides a simple model of decoherence which exhibits both relaxation 
and dephasing. Focusing on the semiclassical limit, we derive equations of motion for the evolution 
of a coherent state and investigate the evolution of a superposition of such states. While a coherent 
state evolves on a timescale given by the classical ratio of the angular momentum of the gyroscope 
to that of the incoming particles, dephasing acts on a much shorter timescale that depends only on 
the angular difference of the states in the superposition. 



I. INTRODUCTION 

A quantum gyroscope is a quantum system that singles 
out a direction in space. It allows for the measurement of 
the spin of incoming particles through the measurement 
of the total spin of both the gyroscope and the incom- 
ing particle. This is one important example of the more 
general theory of quantum reference frames [J. While 
the gyroscope allows the measurement of spin along its 
direction, it will necessarily degrade with time. We will 
briefly review work that has been done in that direction 
before revisiting the same model but not as a quantum 
reference frame, but as a simple model of decoherence. 

In 0, Bartlett, et al. studied a model in which the 
reference frame is represented by a spin € in a polarized 
state. A succession of unpolarized spin- 1/2 particles (all 
initially in the same state) interact with the reference, the 
interaction chosen being a measurement of the total spin, 
after which the particle is discarded and the result of the 
measurement is ignored. The number of interactions acts 
as a time parameter, and they showed that this "quantum 
gyroscope" would degrade on a time scale proportional to 
P . A similar model was studied by Poulin and Yard 
who considered the spin- 1/2 particles to be partially po- 
larized along a fixed axis. They found that, in addition 
to degradation, the reference spin drifted towards the di- 
rection of the polarization of the spin- 1/2 particles on a 
time scale linear in £, making it in a sense a more danger- 
ous effect than degradation since, for large £, it occurs on 
a faster time scale. They concluded that the gyroscope 
would evolve semiclassically and be useful to measure 
spin-i particles along the drifted direction rather than 
its original polarization. In Ahmadi, et al. analyzed 
the evolution of the gyroscope after repeated measure- 
ments and showed that by retaining the outcomes of the 
measurements, one could correct the drift of the quantum 
gyroscope. 

In this paper, building on previous work we use 
a similar model of quantum gyroscope but focus on its 



decoherence in order to analyze the transition between 
quantum and semiclassical behavior. From this new 
point of view, the model is promising since it exhibits 
a sharp transition between these two regimes. The setup 
we study is representative of many decoherence processes: 
a quantum system placed in a semiclassical state (the gy- 
roscope) interacts with a large number of incoming parti- 
cles which together form its environment. However, un- 
like most decoherence models, the transfer of information 
to the environment occurs since each incoming particle 
interacts successively with the quantum gyroscope. 

Another important motivation for this work is linked 
to a strategy recently devised to reduce decoherence in 
quantum dot spin qubits. In these systems, it was real- 
ized that the main source of decoherence is the hyperfine 
interaction between the spin qubit and the Overhauser 
field h resulting from all the nuclear spins in the sub- 
strate. The corresponding Hamiltonian is if = S-h where 
S is the spin operator of the spin qubit. In the presence 
of an external magnetic field, the essential contribution 
to the interaction Hamiltonian is due to the component 
of the spin along the direction of the magnetic field. The 
interaction thus reduces to a tensor product. For such a 
Hamiltonian preparing an eigenstate of the Overhauser 
field h will greatly reduce decoherence 043 j and will 
even suppress it if the prepared state is also an eigen- 
state of the self-evolution of the nuclear spins . There- 
fore, recent experiments have tried to pola rize the nuclear 
bath [lol . [Tl| or lock its polarization |l2l . [l3| . Preparing 
the initial state of the nuclei in a fully polarized state 
would be exactly the situation we consider in this paper 
if the nuclear spins were restricted to be in a symmetric 
state and could therefore be considered as one particle 
with large spin Therefore, the model we use can be 
considered a toy-model that describes the relaxation and 
dephasing of the nuclear spins when they interact with 
stray electrons. 

In the next section, we describe the model studied 
which analyzes the effect on the quantum gyroscope of 
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a succession of interactions with spin-1/2 particles. This 
effect can be written as a quantum channel which, with a 
particular choice of the interaction time, is very similar to 
the quantum channel where a joint measurement of the 
total angular momentum is performed In Section III 
we obtain equations for the orientation and magnitude 
of the average value of the angular momentum of the 
gyroscope in the semiclassical limit, and discuss quan- 
tum effects. In Section IV we discuss the evolution of 
a coherent state and a superposition of two such states. 
In the fifth section we argue that coherent states should 
minimize purity loss but nonetheless suffer a significant 
purity loss initially. In the final section, we discuss what 
can be learned about the state of the reference by mea- 
suring the spin of the interacting particles rather than 
simply discarding them. 



II. THE MODEL 

Our model is summed up in the following timeline. At 
t = 0, a, spin-^ particle S in a state described by a den- 
sity matrix ^ starts to interact with the spin-^ gyroscope 
TZ according to a Heisenberg Hamiltonian 'Hns = L ■ S 
where L (resp. S) is the angular momentum operator of 
the gyroscope (resp. of the particle). At t ~ t, the inter- 
action ceases and the particle is discarded. This defines 
a quantum channel £ mapping the state of the gyroscope 
TZ before interacting with the particle, represented by a 
density matrix p of dimension d — 21+1, to its state after 
the interaction. Immediately thereafter, a second parti- 
cle, also in state ^, starts to interact with the gyroscope, 
the Hamiltonian and duration as above. This process is 
then repeated n times. 

This model has also been discussed by Ahmadi, Jen- 
nings and Rudolph recently Q; their main goal was to 
see if the drift of the gyroscope frame noted in Q could 
be eliminated. They considered and compared several 
possible strategies to correct for this drift. 

We can write a discrete iterative equation for the state 
of the gyroscope. If after n interactions it is in the state 
(where £°{p) = p), then 



£"+^(p) ^ p((n + l)r) 



= Trs [ 



(1) 



The evolution operator e"*^'^'^'^ can be expressed an- 
alytically by noting that V fc, (L • S)'' = Ofel -I- 6fc(L • S). 
The coefficients and 6^ obey a set of coupled linear 
recurrence relations 



Ofc+i — -, Ok 



Ofe+1 = ak - —■ 



(3) 
(4) 



a(T), &(t) are given by 
air) 



d + 1 _jd-l d~l jd+l 

e + —T^e'—^, (5) 
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b{r) = t (e-^^ - e'^^) 



(6) 



This technique, using a decomposition of the k*^ power 
of the interaction Hamiltonian in a finite sum of opera- 
tors, might be used for interaction with spins of higher 
dimension. 

Finally, the calculation of ([2]) for oneparticle interact- 
ing with the gyroscope {n ~ 0) yields [2, Q 
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Thus, e"*^^'^'^ = a(T)I -f- 6(t)L ■ S, where the functions 



This quantum channel, for interaction time t = , yields 
a similar quantum channel than the one based on measur- 
ing the total angular momentum described in Q, except 
for the last term (absent in Q), which corresponds to a 
precession around the axis of preferred polarization (S). 
Thus, up to a slight rotation, interaction and joint mea- 
surement imply the same evolution of the gyroscope. It 
thus seems that this interaction time maximize the ex- 
traction of information on the incoming particles and we 
thus expect it to maximally degrade the gyroscope. 

In order to pursue the analysis further, we will choose 
coordinates as in Q such that the z axis corresponds 
to the polarization of the incoming particles and the x 
axis is chosen so that initially (L) lies in the xz plane. 
We assume ^ = ^ + 2{Sz)Sz. Rather than working with 
operators L it is convenient to rotate the axes by an angle 
d{t) around the y axis such that = cosd — s'md Lz 
and = s'mO Lx + cos6 Lz- The angle 9{t), represented 
in Fig. [2 is chosen so that {L^'-^Ht)) = and {Ll'-^\t)) 
ir{t), where < r{t) < 1. 



III. SEMICLASSICAL EQUATIONS OF 
MOTION 

A. Systematic method 

Equations of evolution for 9{t) and r{t) would allow 
us to determine the quasiclassical evolution of the gyro- 
scope. An equation for 9{t) was already derived in Q, 
but not a corresponding equation for r{t). Here, we pro- 
vide a systematic method that allows us to derive an 
equation for r{t) which (like that for d{t)) is valid to or- 
der 0(l/£2). More specifically, we show that macroscopic 
expectation values (L^) and (I/^) at time t + 1 can be ex- 
pressed analytically in terms of ((L^) }, ((L^) }, ((i^) ) 
and (L^Ly) at time t, where 9 is the angle corresponding 
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Figure 1: The axes are defined sucii that the gyroscope points 
in the Q{t') direction witli respect to the z axis and the po- 
larization of the gyroscope is £r{t). 6{t) naturally defines a 
rotated system of axes, labeled xg and zg. 

to the choice of axes at time t. To do so, we write the 
quantum channel ([7]) as a function of the by using the 
Kraus form provided in Q. 

+ P{S.)sme{LlpLl~LlpLl) 
+ ^z{S.^ cose {LlpLl-LlpLl) 
+ 1{S,} cose (Llp + pLl) 
-1{S,) sine {Lip + pLi) (8) 
Calculation of {Lf{t + 1)) = Tr [Lf yields 

{Ll{t+1)) = -±{S.)cose{LlLl) 

~^{S.)sine[iii(+l)-{{LlY)){9) 
{Ll{t + l)) = il-^){Ll) + ±{S,)sme{LlLl} 

+ ^{S^)cose[iie+l)-{{Ltf))ilO) 

These expressions will allow us to derive equations of 
motion in the semiclassical regime. 



regime, we assume 

Llp^erp (11) 

which implies that ((L^)^) ~ i'^r^ and (L^L^) ~ 0. 
Furthermore, we are interested in the semiclassical limit 
£ ^ I. In that limit, numerical simulations show that 
~ {{LlY) ~ {e{£+l)-£\^) 12. The average 
values thus reduce to 

(L^(t + 1)) ^ -i±!f(5,)sin0 (12) 
(Lf(t+1)) ~ fr + (l-r2)(5^)cos0 (13) 

We can then derive the equations of motion on the 
macroscopic variables r(t), 0(t) by using 

e ^ (i«,(t + i))/(Lf(t + i)) (14) 

2rf = r^{t + \)-r'-{t) (15) 

The final result is 

f = A (l-r2(i))cos6'(i)+ 0(1/^2) (16) 

^ = -A^^sin^?(t)+0(l/£2) (17) 

where A — {Sz)/£ sets the time scale of the evolution. 

Equations ([T6|) and (|17p provide a clear picture of the 
semiclassical evolution of the gyroscope. This set of equa- 
tions improves on the work in Q which only derived 
heuristically a differential equation for e{t) valid in the 
limit r = 1. Our method is more systematic and allows 
us to derive an equation for r{t). Furthermore, in it 
was argued that r remains close to 1 based on numeri- 
cal observations. We strengthen this statement since Eq. 
()17p provides an analytic confirmation that r does indeed 
remain close to 1 when r(0) ~ 1 and 7^ tt. 

The striking feature of the set of coupled equations 
(fTHll and (fT7|) is that the time scale is entirely set by the 
ratio A of the angular momentum of the incoming par- 
ticle {Sz) and the angular momentum of the gyroscope 
£. Thus, the semiclassical evolution of the quantum gy- 
roscope only depends on three parameters: £ and r that 
characterize the semiclassical state of the gyroscope and 
the ratio A. This set of parameters is extremely small 
when compared to the coefBcients necessary to express 
the states of the gyroscope and the incoming particle as 
vectors in their respective Hilbert spaces. From an op- 
erational point of view the semiclassical evolution is en- 
tirely described by these three parameters and does not 
depend on the microscopic details of the gyroscope and 
the incoming particle. 



Evolution from a semiclassical state 



C. Quantum effects for a finite gyroscope 



Now, suppose that p is in a semiclassical state, i.e., a 
state that is close enough to an eigenstate of . In this 



Numerically, one notices that a gyroscope initially 
fully polarized will undergo an initial polarization drop 
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as can be seen on Fig. [31 This phenomenon might 
shed new light on recent experiments involving quantum 
dots [lol . [ill where the nuclear bath reached a level of 
polarization lower than had been expected. The semi- 
classical equation of motion ([TB)) alone fails to explain 
why the polarization undergoes an initial drop. Indeed, 
for a quantum gyroscope initially fully polarized (r = 1), 
eq. (fTB)) breaks down because the evolution of the po- 
larization is dominated by corrections to the otherwise- 
dominant term. A more careful derivation yields 



xlO" 



Afl-r2)cos6'- 



1 + {S.)' A 



2^2 



— sin^ Or-^ (18) 



in the limit where i ^ 1 and for 9 > it/ 2. When the 
gyroscope is fully polarized and 9 is close to tt, the sec- 
ond term of (jl8p dominates and the polarization decays 
exponentially which explains the initial drop that brings 
the polarization to the form 1 — £. The first term of ([18]) 
takes over but vanishes again for 6 = tt/2. Thus cor- 
rections are responsible for both the initial drop and the 
precise position of the minimum of the polarization. 

For small 9 eq p6p is adequate and its linearization 
yields 



-2Acos 9£. 



(19) 



Assuming that the loss of polarization remains small 
enough, we can use the semiclassical (r = 1) solution Q 
of ([IZD 



tan [9{t)/2] 



tan [9{0)/2] 



(20) 



that shows that the gyroscope aligns itself with the axis 
of polarization of the incoming particle at a rate given by 
the semiclassical ratio A. Injecting this in the linearized 
equation of motion, we can integrate cos by expressing 
it in terms of tan [9{t)/2] and using ([20)) to turn it into a 
rational expression of e~^*. The result reads 



e{t) = e(T) 



(cosh At + cos 9q sinh At) 
(cosh Xt + cos 9o sinh At)^ 



(21) 



We compare these analytical expressions with numer- 
ical simulations. Numerical results are summarized in 
Fig. [21 four regimes of behavior appear, delineated by 
times ti, <2 and t^. For t < ti, corrections of higher or- 
der than 1 play an important role and Eq. ([TS]) is not 
exact, as can be seen in the top figure. Thus, the initial 
drop is not captured by semiclassical equations. How- 
ever, Eq. ([T5)) accounts well for the polarization decrease 
for ti < t < t2. However, while a formal solution can 
be written down by linearizing the equation and using 
the method of variation of parameters, we failed to at- 
tain a compact analytical formula. We thus performed a 
numerical simulation of Eq. (|18|) which agrees well with 
the brute force computation of r from the quantum state 
evolving under the quantum channel. For the crossover 
period t2 < t < t^, neither (|18p nor ([TB)) accounts for the 
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Figure 2: Comparison between the evolution of the polariza- 
tion r{t) by brute force numerical simulation of the quantum 
channel and the analytical solution. Four regions appear. The 
initial drop [t < ti) and the crossover region {t2 < t < t^) are 
not described by the semiclassical eq. of motions and seems 
to require a full quantum treatment. On the contrary, the 
polarization decrease {ti < t < t2) and increase (ts < t) are 
described very accurately by semiclassical equations. 



evolution of r. A better assumption than (|lip is prob- 
ably needed to properly describe the transition between 
([T5)) and ([TB]) . The return to full polarization for t > 
is well described by the analytical expression ([?!]) with 
initial condition t = t^. 



IV. DECOHERENCE OF A SUPERPOSITION 
OF COHERENT STATES 

In this section, we will describe the evolution of states 
that behave semiclassically, namely coherent states, be- 
fore investigating the evolution of a superposition of such 
states. It will exhibit dephasing on a much shorter 
timescale than relaxation. 



A. Coherent states 

Coherent states have maximal angular momentum 
along a certain axis. If this axis is in the xz plane with 
polar angle 9, the corresponding quantum state reads, in 



the |£, m) basis, 



{9\£m{om 



\£,£). 



(22) 



The semiclassical equations of motion (|T6|) indicate that 
coherent states will keep maximal polarization (Vt r{t) ^ 
1) and align themselves with the z axis, rotating around 
the y axis according to (|20p. Once again, the ratio A 
plays a crucial role as it sets the timescale of the align- 
ment of the gyroscope with the z axis, which is expo- 
nentially fast. Equation ((20| describes the relaxation or 
thermalization of the gyroscope, a classical phenomenon 
governed by the semiclassical equation of motions. In 
the next paragraph, we will investigate a purely quan- 
tum phenomenon, dephasing, which typically acts on a 
much shorter timescale. 



B. Evolution of a coherent superposition of 
coherent states 

Consider a gyroscope initially prepared in a superposi- 
tion of coherent states jV') = a\9) + b\cf>). The evolution of 
the coherence terms can be computed by usingthe Kraus 
expression of the quantum channel given in [3| and ex- 
pressing operators acting on \9) (resp. in the basis 
{Ll, Ll) (resp. {L^, Lf}). The following equation de- 
scribes how the evolution reduces the coherence terms: 



(|0)((/)|)|0)=cos2. 



COS 6 



(cos ( 



2{S,) - cos^ 



(23) 



In the semiclassical limit, i.e.. 1 the first term dom- 
inates, as confirmed by numerical simulations shown in 
Fig. [3]for (/i = 0. 

Thus, the coherence terms vanish more and more 
rapidly as the initial value of \6 — 0| increases. This is a 
feature similar to the result obtained for Quantum brow- 
nian motion [l5| where the coherence time of a super- 
position of two localized wavepackets is inversely propor- 
tional to the square of the distance [l^. The difference 
of angles characterizes the classical distance between the 
two quasi-classical states. 

The superposition thus evolves according to 



£{\i)){il^\) = \a\''\e + de){e + de\ + \b\''\(j)- 



+ ah* cos 



|(?)(0|-fX +h-c. (24) 



While eq. (|23|) gives the expression of a matrix ele- 
ment, eq. is stronger since it represents a density 
matrix. However, we had to introduce an operator x, 
required by the trace-preserving nature of the channel £. 
This operator accounts for terms that do not play a role 
in subsequent application of the channel evolution in the 
semiclassical limit ^ ^ 1 (see Fig. [5]) as pointed out by 
numerical evidence shown on Fig. 21 




Figure 3: Evolution of the off-diagonal term \0){Q\ correspond- 
ing to the superposition a\9) -f /3|0) where |0) = points 
in the direction z and \9) is given by equation (|22|l 
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Figure 4: Numerical evolution of the trace norm of the oper- 
ator X with different values of £ for 6 = n/?) (top) and with 
different values of 9 for ^ = 20 (bottom). The norm vanishes, 
indicating that x does not play a role in the evolution of the 
superposition of coherent states. 



Thus, IV') (V" I turns into a mixture 



6p|(/)(t))(0(t)| after interacting with a few incoming par- 
ticles. This decoherence process is due to the entangle- 
ment with the incoming particles that allows the envi- 
ronment (made of all the incoming particles) to partially 
distinguish the two coherent components {\0) and \(f))) of 
the state. In the semiclassical limit d ^ 1, decoherence 
acts on a timescale that essentially depends on the "clas- 
sical" distance — 0| and not on the semiclassical ratio A 
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Figure 5: Evolution of the the purity Trp^{t) for increasing 
values of £. Dashed lines represent the loss of polarization 
e — 1 — r and the solid lines represent the inverse of dTip^. 
Before the loss of polarization reaches a maximum, solid and 
dashed lines agree for all values of £. 



initially. After a few interactions, the gyroscope will lose 
a small amount of polarization, i.e. r{t) = 1 — e where 
e is small for large I. The resulting state p will then be 
roughly a statistical mixture of ed states among the d 
available states. Supposing that each coefficient has the 
same order of magnitude, the purity Trp^ will scale as 
£d{-^Y = ^7d^ explaining the loss of purity. This predic- 
tion agrees with numerical results (Fig. [5]). 

VI. INDUCED POVM 

Interaction between the quantum gyroscope TZ and the 
spin-i particle, which can be represented by a unitary op- 
erator U = e~*^^'5'^, creates an entangled state. There- 
fore, a measurement represented by a POVM {n±} on 
the particle yields information about the state of the gy- 
roscope according to the POVM {A±} 

A± = Trs [(Ik ViWHln «> n±)f/(% V?)] (25) 

Our calculations show that a projective measurement of 
the spin-i particle along an axis u yields 



that governs relaxation. Indeed, unless the two coherent 
components are very close (6* ~ (/)), decoherence acts ex- 
ponentially fast, in the sense that {9\£ {\0){4>\) \<j>) reduces 
to a value e after interacting with a number of incom- 
ing particles that scales as — log ^cos^ ^T^) ■ Thus, in 
general, decoherence occurs on a timescale much shorter 
than relaxation. Once the two components \9) and |0) 
have decohered, they each relax and align with the z axis 
independantly. These features are a revealing example of 
characteristics expected from the general theory of deco- 
herence. 



V. LOSS OF PURITY 

The predictability sieve, introduced in identifies 
semiclassical states as those which minimize entropy pro- 
duction or purity loss. In the same way that minimal 
uncertainty coherent states are semiclassical for quantum 
Brownian motion m our model, the SU{2) coherent 
states clearly behave semiclassically and should minimize 
purity loss. 

Interestingly, they still undergo a significant purity loss 



A± = ^Tts [n±C] In ± ((S) X u) • ^ (26) 

corresponding to a measurement of L along (S) x u. 

A broader problem is to evaluate the POVM induced 
by n spin-i particles interacting one after the other with 
the gyroscope during a time r and then measured col- 
lectively. It is not clear how much information about 
the quantum state p of the gyroscope can be gained this 
way. It is possible that only information about a few 
global properties of p can be extracted for large but fi- 
nite n thus making these properties good candidates for 
characteristics that emerge in the semiclassical limit and 
which correspond to macroscopic attributes of the gyro- 
scope. 
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